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In  the  present  study,  the  non-linear  solar  radiation  models  for  predicting  the  monthly  average  daily  diffuse 
radiation  are  developed  using  the  measured  data  on  global  radiation,  diffuse  radiation  and  sunshine  hours 
for  12  locations  of  India.  Statistical  method  is  used  to  derive  these  correlations.  The  developed  models  are 
employed  to  estimate  the  monthly  average  daily  diffuse  radiation.  The  performance  of  these  correlations 
is  compared  with  existing  model.  Accuracy  of  developed  relationships  is  also  tested  using  statistical 
indicators  viz.  Percentage  error  (PE),  root  mean  square  error  (RMSE),  mean  percentage  error  (MPE)  and 
mean  bias  error  (MBE).  The  study  finds  that  these  statistical  parameters  have  very  low  values  for  the 
proposed  models.  A  cubic  correlation  of  diffuse  coefficient  with  percent  possible  sunshine  gives  the  best 
fit.  The  maximum  values  of  RMSE,  MPE  and  MBE  for  the  proposed  third  order  equation  are  4.33%,  8.68% 
and  -1.25%  respectively  while  in  the  case  of  existing  model  these  values  are  13.28%,  13.39%  and  -3.83% 
respectively.  Hence,  it  is  possible  to  apply  the  cubic  equation  for  the  prediction  of  monthly  mean  daily 
diffuse  radiation. 
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1.  Introduction 

The  power  obtained  from  solar  radiation  reaching  the  Earth  is 
many  times  greater  than  the  power  generated  by  man  by  other  con¬ 
ventional  sources.  If  it  is  trapped  in  a  useful  and  effective  manner, 
all  the  present  and  future  human  energy  needs  may  be  met  on  a 
continual  basis. 

The  energy  received  from  the  Sun  is  absorbed  and  also  reflected 
back,  a  part  of  which  is  trapped  in  the  atmosphere  of  the  Earth. 
At  any  place  the  sum  of  direct  and  diffuse  radiation  is  known  as 
global  or  total  radiation.  Solar  energy  can  be  utilized  for  different 
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applications  by  a  suitable  choice  of  tapping  device.  The  availability 
of  solar  radiation  will  depend  on  place,  on  tilt  of  the  surface  on 
which  it  is  received  and  also  upon  a  number  of  other  factors. 

Since  the  solar  radiation  reaching  the  Earth’s  surface  depends 
upon  climatic  conditions  of  the  place,  a  study  of  solar  radiation 
under  local  climatic  conditions  is  necessary.  This  knowledge  of 
local  solar-radiation  is  essential  for  the  proper  design  of  buildings, 
energy  systems,  solar  energy  systems  and  a  good  evaluation  of  ther¬ 
mal  environment  within  buildings  [1-6].  It  is  measured  in  terms  of 
sunshine  duration  or  in  terms  of  direct,  diffuse  and  global  radiation. 
In  the  absence  of  solar  radiation  data,  the  need  for  an  empirical 
model  for  prediction  of  solar  radiation  is  desirable. 

Solar  radiation  models  are  of  two  categories  namely  paramet¬ 
ric  model  in  which  detailed  information  of  atmospheric  conditions 
is  required  and  decomposition  models  which  require  global  radi¬ 
ation  for  the  estimation  of  direct  and  diffuse  component  [7].  The 
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metrological  parameters  such  as  amount  and  distribution  of  clouds, 
fractional  sunshine,  atmospheric  turbidity  and  perceptible  water 
content  are  used  in  parametric  models  [8]. 

Several  empirical  correlations  have  been  developed  to  estimate 
solar  radiation  [9-16].  Orgill  and  Hollands  [17]  have  used  the  fol¬ 
lowing  diffuse  radiation  model  based  on  the  hourly  data  of  Toronto 
(Canada)  correlating  the  diffuse  fraction  Kd  and  clearness  index  KT: 

77  =  1  -  0.249 /<r,  Kt  <  0.35 
Hg 

77  =  1.577  -1.84 Kt,  0.35  <  KT  <  0.75 
Hg 

77  =  0.177,  Kj  >  0.75  (a) 

Hg 

Erbs  et  al.  [18]  studied  the  data  from  location  of  USA  and  have 
given  the  relations  as: 

77  =  1-0.09 Kt,  Kj  <  0.22 
Hg 


77  =0.61092  +  3.6259  KT  -  10.171  K?  +  6.338  K?, 

Hg 

0.2  <  /CT  <  0.7 

77  =  0.672  -  0.474 Kt,  KT  >  0.7  (f) 

Hg 

Studies  [22-24]  have  also  been  made  to  estimate  hourly  diffuse 
solar  radiation  based  on  artificial  neural  network  technique.  Liu 
and  Jorden  [25]  have  developed  model  for  estimating  the  diffuse 
components  of  daily  global  radiation  which  is  given  by: 

77  =0.384  -0.41  6KT,  (g) 

Spencer  [26]  studied  the  mean  daily  data  from  the  localities 
in  Australia  and  proposed  following  latitude  (0)  dependent  rela¬ 
tionship  between  diffuse  fraction  of  global  radiation  and  clearness 
index. 

77=a-b  Kj,  0.35  <  I(j  <  0.75 
Hg 


77  =0.9511  -0.1604KT  +  4.388  K^  +  16.638  K|  +  12.336  id, 
Hg 

0.22  <  Kj0.22  <  0.8 


77  =  0.165,  Kt>  0.8  (b) 

Reindl  et  al.  [19]  proposed  two  models  for  hourly  diffuse  radia¬ 
tion  using  the  data  of  USA  and  Europe.  The  first  model  is  in  terms 
of  clearness  index  and  is  given  by: 

77  =  1-02  -  0.248  Kt,  Kt  <  0.3 
Hg 

77  =  1 .45  -  1 .67  Kj,  0.3  <  Kj  <  0.78  (c) 

Hg 

The  second  model  has  been  given  in  terms  of  clearness  index 
and  solar  elevation  a  and  is  as  follow: 

^7  =  1. 02  -  0.248 /CT  + 0.01 23  sin  a,  KT  <  0.3 
Hg 

77  =  1.4-  1.749 Kt  +  0.177  sin  a,  0.3  <  Kj  <  0.78 

77  =  0.486 Kt-  0.182  sin  a,  KT  >  0.78  (d) 

™g 

Lam  and  Li  [20]  have  developed  following  relationship  for  the 
prediction  of  hourly  diffuse  components  of  global  radiation  for 
Hong  Kong. 

77  =  0.977,  Kt  <  0.15 

™g 

77  =  1.237-  1.361  JCT,  0.15  <  Kj  <  0.7 

™g 


where, 

a  =  0.94  +  O.O118|0|,  b  =  1.185  +  O.O135|0|,  (h) 

Boland  et  al.  [27]  have  given  a  logistic  function  of  the  form 
/(x)  =  a0/l  +a^ea2x  for  estimating  diffuse  solar  radiation  for  Aus¬ 
tralian  condition  and  Jacovides  et  al.  [28]  have  proved  its  validity 
for  locations  in  Cyprus. 

The  objective  of  the  work  reported  in  this  paper  is  the  develop¬ 
ment  of  non-linear  relations  for  estimating  the  diffuse  component 
of  global  radiation  on  a  horizontal  surface  and  to  compare  these 
relationships  with  the  calculations  based  on  existing  model  and 
then  to  find  out  the  most  suitable  correlation  for  computing  diffuse 
radiation  for  locations  in  India.  The  monthly  average  data  of  solar 
radiation  are  employed  in  partial  regression  analysis  to  derive  the 
models. 


2.  Methodology 


Despite  using  different  approaches  almost  all  the  authors  viz. 
Page  [29],  Gupta  et  al.  [30],  Modi  and  Sukhatme  [31],  Mani  and 
Rangrajan  [32]  and  Gopinathan  [33]  have  given  following  linear 
equation  for  the  estimation  of  monthly  mean  daily  diffuse  radia¬ 
tion: 


Hd 

Hg 


=  a  +  b 


Hg\ 

Ho) 


where  Hd  is  the  monthly  mean  daily  diffuse  radiation  on  a  hor¬ 
izontal  surface  in  kWh/m2,  Hg  is  the  monthly  mean  daily  global 
radiation  on  a  horizontal  surface  in  kWh/m2  and  H0  is  the  monthly 
mean  of  extra-terrestrial  radiation  in  kWh/m2.  a  and  b  are  correla¬ 
tion  coefficients. 

In  our  study,  the  experimental  data  on  global  radiation,  diffuse 
radiation  and  sunshine  hours  reported  by  Mani  and  Rangarajan  [34] 
have  been  analyzed  for  14  locations  of  India.  We  have  taken  a  group 
of  following  relation: 


77  =  0.273,  Kt  >  0.7  (e) 

Hg 

Mondel  et  al.  [21  ]  have  presented  diffuse-global  model  using  the 
hourly  data  for  Northen  Ireland  as: 


(i)  models  correlating  diffuse  fraction  Hd/Hg  with  clearness  index 
I<j  and  fraction  of  possible  number  of  sunshine  S/Smax 


Hd 

Hg 


a  +  b 


77  =  1.237-  1.361  Kx, 
Hg 


Hd 

Hg 


a  +  b 


+  c 


Hg\ 

H0) 


2 


(3) 
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Table  1 

Regression  coefficients  for  the  chosen  locations. 


Stations 

For  third  order  equation 

For  second  order  Equation 

a 

b 

c 

d 

a 

b 

c 

Ahmadabad 

0.5456 

-1.1613 

1.744 

-0.9945 

— 

— 

— 

Bhavnagar 

0.308 

0.068 

-0.078 

-0.155 

0.28 

0.236 

-0.37 

Mumbai 

0.3509 

-0.4695 

1.0988 

-0.9116 

- 

- 

- 

Kolkata 

0.3401 

-0.3815 

0.8308 

-0.6857 

0.2378 

0.2387 

-0.3323 

Goa 

0.2537 

0.1737 

-0.0495 

-0.2947 

1.0654 

-0.345 

-0.6604 

Jodhpur 

-0.0684 

1.8254 

-2.7068 

1.085 

0.35 

0.0459 

-0.2672 

Madras 

0.2753 

0.3472 

-0.9128 

0.4683 

0.4195 

-0.3099 

0.0597 

Nagpur 

0.3181 

-0.0986 

0.0422 

-0.1253 

0.2957 

0.0334 

-0.1895 

New  Delhi 

0.7863 

-2.0503 

2.7041 

-1.3029 

0.3887 

-0.2137 

-0.0174 

Pune 

0.1982 

0.8814 

-1.7625 

0.8095 

0.4113 

-0.2313 

-0.0575 

Tiruvanantpuram 

0.3751 

-0.2126 

0.0570 

-0.0980 

0.3350 

-0.0344 

0.1838 

Vishahapatanam 

0.3503 

-0.1720 

-0.0068 

-0.0119 

0.3426 

-0.1405 

-0.0435 

Hd 

Hg 


a  +  b 


+  c 


Hg\ 

Ho) 


+  d 


Hg 

Ho 


3 


Hd 

Hg 


a  +  b 


+  c 


Hg 

Ho 


2 

T  d 


Hg 

Ho 


(4) 

(5) 


The  accuracy  of  these  relationships  is  also  tested  by  calculating 
RMSE  and  MBE.  These  are  defined  as: 


RMSE  = 


E(Hi, 


Hi,0)2] 


1/2 


n 


(11) 


Hd 

Hg 


a  +  b 


(ii)  models  correlating  diffuse  coefficient  Hd /H0  with  fraction  of 
possible  number  of  sunshine  S/Smax 


Hd 

Ho 


=  a  +  b 


+  c 


2 


Hd 

Ho 


=  a  +  b 


(7) 

(8) 


mbe=[£(H!£ — (12) 

n 

where  H[  c  is  the  ith  calculated  value,  H[  o  is  the  ith  measured  value 
and  n  is  the  number  of  observations. 

RMSE  provides  information  on  the  short-term  performance  of 
a  model.  The  lower  the  RMSE,  the  more  accurate  is  the  estimation 
[21  ].  MBE  provides  information  on  the  long-term  performance  of  a 
model.  Positive  and  negative  MBE  values  show  overestimation  and 
underestimation  respectively  [21]. 

To  check  the  linearity  between  the  estimated  and  observed  val¬ 
ues,  coefficient  of  determination  R 2  is  used. 


Here  S  is  monthly  mean  of  sun-shine  hours  per  day  and  Smax  is 
monthly  mean  of  maximum  possible  sun-shine  hours  per  day  and 
a,  b ,  c,  d,  e  and/are  the  regression  coefficients.  The  no.  of  regression 
coefficients  changes  with  the  order  of  the  equation. 

We  have  used  partial  regression  analysis.  The  experimental  data 
of  each  station  are  used  in  normal  equations  and  statistical  analy¬ 
sis  is  carried  out  to  find  the  regression  coefficients  from  the  above 
relations  (from  Eqs.  (2)  to  (8))  for  each  station.  The  regression  coef¬ 
ficients  of  Eq.  (1)  are  also  computed  for  each  of  the  12  locations.  A 
comparison  is  then  made  between  the  results  obtained  on  the  basis 
of  Eq.  (1)  and  above  proposed  models  (from  Eqs.  (2)  to  (8)). 

3.  Statistical  parameters 

The  statistical  parameters  have  been  used  to  evaluate  the  per¬ 
formance  of  the  formerly  given  different  models  (Eqs.  (l)-(8))  for 
each  station.  These  are  -  percentage  error  (PE),  mean  percentage 
error  (MPE),  root  mean  square  error  (RMSE)  and  mean  bias  error 
(MBE).  We  have  expressed  RMSE  and  MBE  as  percentage. 

The  percentage  error 

PE  =  ( Hd°~ Hdc )  x  1 00  (9) 

V  Hdo  ) 

is  calculated  for  each  month  of  the  year  for  all  the  proposed  models 
and  also  for  relation  given  by  Eq.  ( 1 ).  Hg0  is  the  observed  and  Hgc  the 
calculated  value  of  diffuse  radiation.  The  mean  percentage  error  is 
also  computed. 

ry(Hdo-Hdc/Hdo)xl00l  , 

MPE  =  ^  -  '  -  (10) 

n 


4.  Results  and  discussions 

Using  partial  regression  analyses,  the  data  on  Hd,Hg  and  S  for  the 
12  locations  are  analyzed  and  statistical  analyses  is  carried  out.  The 
comparison  between  experimental  and  computed  values  of  diffuse 
radiation  from  Eqs.  (2)-(8)  is  made.  We  found  that  for  Ahmad- 
abad  and  Mumbai  model  8  gives  better  results  than  that  of  other 
equations.  So  for  these  two  stations  cubic  relation  between  dif¬ 
fuse  coefficient  ( Hd/H0 )  and  fraction  of  possible  number  of  sunshine 
(5/Smax)  (Eq.  (8))  must  be  used. 

Good  agreement  between  the  observed  and  estimated  diffuse 
radiation  from  Eqs.  (7)  and  (8)  and  is  found  for  the  stations  Bhav- 
nagar,  Kolkata,  Goa,  Jodhpur,  Madras,  Nagpur,  New  Delhi,  Pune, 
Tiruvanantpuram  and  Vishakhapatanam,  in  comparison  to  the 
other  model  (i.e.  Eq.  (1 )).  Hence  square  and  cubic  relation  between 


Table  2 

RMSE  (%)  for  the  proposed  cubic  and  square  correlation. 


Stations 

Third  order 
equation 

Second  order 
equation 

Eq.  (1) 

Ahmadabad 

1.73 

— 

10.4 

Bhavnagar 

0.58 

0.58 

5.2 

Mumbai 

0.58 

- 

8.08 

Kolkata 

4.33 

4.91 

9.53 

Goa 

2.31 

2.31 

4.62 

Jodhpur 

2.6 

2.6 

8.37 

Madras 

0.29 

0.29 

2.6 

Nagpur 

2.31 

2.6 

3.75 

New  Delhi 

2.89 

2.31 

13.28 

Pune 

2.89 

2.89 

7.22 

Tiruvanantpuram 

0.87 

0.87 

2.31 

Vishahapatanam 

0.87 

0.87 

1.15 

/.  Karakoti  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  15  (201 1 )  2378-2384 


2381 


Fig.  1.  Comparative  plot  of  the  percentage  error  for  the  proposed  models  and  for  Eq.  (1)  for  Ahmadabad,  Bhavnagar,  Mumbai,  Kolkata,  Goa  and  Jodhpur. 


Fig.  2.  Comparative  plot  of  the  percentage  error  for  the  proposed  models  and  for  Eq.  (1 )  for  Madras,  Nagpur,  New  Delhi,  Pune,  Tiruvanantpuram  and  Vishakhapatanam. 
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Fig.  3.  RMSE  (%)  for  the  proposed  third  order  equation,  second  order  equation  and 
for  Eq.  (1). 

diffuse  coefficient  and  percent  possible  sunshine  (i.e.  Eqs.  (7)  and 
(8))  must  be  applied  for  these  locations.  The  regression  coefficients 
of  the  proposed  equations  for  each  station  are  found  out  which  are 
listed  in  Table  1. 

It  is  clear  that  second  and  third-order  equation  between  diffuse 
coefficient  Hd/H0  and  percent  possible  sunshine  S/Smax  are  fit  for 
all  the  considered  stations.  Thus  these  square  and  cubic  relations 
can  be  used  for  majority  of  locations  in  India. 

To  check  the  accuracy  of  these  relationships  plots  of  percentage 
errors  (PE)  for  the  proposed  models  and  that  of  Eq.  (1)  are  drawn 
for  each  location.  These  are  shown  in  Figs.  1  and  2.  For  Ahmadabad, 
Bhavnagar,  Mumbai,  Kolkata,  Goa  and  Jodhpur  percentage  errors 
are  illustrated  in  Fig.  1  while  for  the  other  six  locations  Madras,  Nag¬ 
pur,  New  Delhi,  Pune,  Tiruvanantpuram  and  Vishapatanam  these 
are  shown  in  Fig.  2. 

RMSE  (%)  for  the  proposed  cubic  equation,  second  order  equa¬ 
tion  and  that  of  Eq.  (1)  is  tabulated  in  Table  2  and  plotted  in  Fig.  3. 
Mean  percentage  error  (MPE)  for  these  two  models  along  with  the 
widely  used  relationship  (Eq.  (1 ))  is  illustrated  in  Table  3  and  Fig.  4. 
Plots  are  also  drawn  for  percentage  MBE  (listed  in  Table  4)  which 
are  shown  in  Fig.  5. 

For  the  12  locations,  the  percentage  errors  as  shown  in 
Figs.  1  and  2  give  the  comparison  between  PE  from  the  proposed 
models  and  from  the  relation  given  by  Eq.  (1).  Fig.  1  indicates  that 
percentage  error  for  calculations  based  on  Eq.  (1)  are  higher  com- 


Table  3 

MPE  for  the  proposed  cubic  and  square  correlation. 


Stations 

Third  order 
equation 

Second  order 
equation 

Eq.(l) 

Ahmadabad 

6.53 

— 

11.5 

Bhavnagar 

0.69 

0.14 

7.1 

Mumbai 

4.45 

- 

9.88 

Kolkata 

5.05 

5.13 

7.97 

Goa 

6.08 

6.05 

7.89 

Jodhpur 

8.68 

8.68 

13.39 

Madras 

4.73 

4.73 

5.54 

Nagpur 

4.45 

4.45 

4.83 

New  Delhi 

5.69 

5.69 

10.53 

Pune 

4.74 

5.2 

10.29 

Tiruvanantpuram 

2.77 

2.62 

5.32 

Vishahapatanam 

5.07 

4.27 

5.83 

Stations 

Fig.  4.  MPE  for  the  proposed  third  order  equation,  second  order  equation  and  for 
Eq.(l). 


Table  4 

MBE  (%)  for  the  proposed  cubic  and  square  correlation. 


Stations 

Third  order 
equation 

Second  order 
equation 

Eq.(l) 

Ahmadabad 

-0.5 

— 

-3 

Bhavnagar 

0.17 

0.17 

-1.5 

Mumbai 

-0.17 

- 

-2.3 

Kolkata 

-1.25 

-1.42 

-2.75 

Goa 

-0.67 

-0.67 

-1.33 

Jodhpur 

-0.75 

-0.75 

-2.42 

Madras 

0.08 

0.08 

0.75 

Nagpur 

-0.67 

-0.75 

-1.08 

New  Delhi 

-0.83 

-0.83 

-3.83 

Pune 

-0.83 

-0.83 

-2.08 

Tiruvanantpuram 

-0.25 

-0.25 

-0.67 

Vishahapatanam 

-0.25 

-0.25 

-0.33 

Proposed  third  order  equation 


Stations 


Fig.  5.  MBE  (%)  for  the  proposed  third  order  equation,  second  order  equation  and 
for  Eq.  (1). 
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Fig.  6.  Coefficient  of  determination  R2  for  the  considered  locations. 


pared  to  the  proposed  models  for  all  the  six  stations  (Ahmadabad, 
Bhavnagar,  Mumbai,  Kolkata,  Goa  and  Jodhpur)  and  in  the  case  of 
Bhavnagar,  these  errors  are  comparatively  high.  Fig.  2  (for  Madras, 
Nagpur,  New  Delhi,  Poona,  Tiruvanantpuram  and  Vishapatanam) 
also  gives  the  similar  results.  Although  for  some  of  the  months  of 
the  year,  errors  from  the  proposed  models  are  slightly  greater  than 
the  widely  used  relation  (Eq.  (1)),  but  on  average,  errors  (%)  from 
Eq.  (1 )  are  significantly  higher  than  that  of  our  models  for  all  the  12 
stations. 

The  RMSE  varies  between  0.29%  (Madras)  and  4.33%  (Kolkata) 
for  the  proposed  cubic  model  while  for  Eq.  ( 1 ) ;  it  has  values  between 
1.15%  (Vishakhapatanam)  and  13.28%  (New  Delhi)  as  clear  from 
Fig.  3.  It  is  found  that  root  mean  square  error  for  the  proposed  third 
order  relation  is  lower  than  that  of  Eq.  ( 1 )  for  all  the  locations.  RMSE 
for  the  second  order  proposed  relation  ranges  from  0.29%  (Madras) 
to  4.91%  (Kolkata).  Again  for  widely  used  relation  (Eq.  (1))  it  has 
higher  values  than  that  of  second  order  equation  as  illustrated  in 
Fig.  3.  The  lower  RMSE  values  of  third  and  second  order  proposed 
equation  demonstrate  their  accuracy  compared  to  Eq.  (1). 

Fig.  4  shows  that  MPE  values  have  the  range  0.69%  (Bhavnagar) 
to  8.68%  (Jodhpur)  for  our  cubic  equation  and  0.14%  (Bhavnagar) 
to  8.68%  (Jodhpur)  for  the  second  order  equation  proposed  by  us. 
While  for  the  relation  given  by  Eq.  ( 1 )  the  values  vary  between  4.83% 
(Nagpur)  and  13.39%  (Jodhpur).  Although  for  some  months  of  the 
year,  percentage  errors  from  proposed  models  slightly  exceed  the 
errors  from  Eq.  ( 1 ),  it  is  clear  from  figure  that  mean  percentage  error 
(MPE)  has  comparatively  lower  values  for  our  proposed  models 
than  that  of  the  widely  used  equation  (Eq.  (1 ))  for  all  stations. 

The  MBE  values  lie  from  0.08%  (Madras)  to  -1.25%  (Kolkata)  for 
the  third  order  equation  and  from  0.17%  (Bhavnagar)  to  -1.42% 
(Kolkata)  for  our  second  order  proposed  equation  as  illustrated  in 
Fig.  5.  For  model  given  by  Eq.  (1 ),  these  values  vary  between  -0.33% 
(Vishakhapatanam)  and  -3.83%  (New  Delhi).  It  is  clear  that  for  Eq. 
(1),  MBE  values  are  higher  compared  to  our  proposed  models  for 
all  stations  which  lead  the  long  term  performance  of  our  proposed 
models. 

The  coefficient  of  determination  R 2  for  the  chosen  locations  is 
plotted  in  Fig.  6.  It  is  clear  from  the  figure  that  R 2  has  a  lower  value 


for  all  the  stations  for  the  proposed  models  compared  to  Eq.  (1). 
The  cubic  correlation  gives  the  best  performance  with  the  highest 
values  of  R2  for  almost  all  the  locations. 

The  statistical  analysis  shows  that  the  best  estimates  of  diffuse 
radiation  are  obtained  from  the  proposed  relationships  compared 
to  the  model  given  by  Eq.  ( 1 ).  For  Bhavnagar,  Kolkata,  Goa,  Jodhpur, 
Madras,  Nagpur,  New  Delhi,  Pune,  Tiruvanantpuram  and  Vishakha¬ 
patanam  both  second  and  third  order  equation  between  diffuse 
coefficient  and  percent  possible  sunshine  give  best  results.  In  case 
of  Ahmadabad  and  Mumbai  only  cubic  equation  gives  best  fit.  Thus 
in  general,  a  third-order  equation  is  found  as  the  most  accurate 
model  for  all  these  locations. 


5.  Conclusions 

The  statistical  parameters  show  that  errors  in  computing  the  dif¬ 
fuse  component  of  global  radiation  from  the  widely  used  relations 
(Eq.  (1 ))  are  higher  than  the  new  proposed  models.  These  proposed 
models  give  good  agreement  between  the  measured  and  computed 
values  of  diffuse  radiation. 

Although  for  most  of  the  locations  in  India,  square  and  cubic 
equations  perform  better  than  the  relation  given  by  Eq.  (1),  third 
order  relationship  between  diffuse  coefficient  and  percent  possible 
sunshine  gives  most  suitable  fit  for  all  the  stations.  Thus  cubic  equa¬ 
tion  can  be  used  to  estimate  the  diffuse  component  of  the  global 
radiation  on  a  horizontal  surface  for  the  locations  inside  India. 

The  models  can  be  applied  for  the  computation  of  local  values  of 
solar  radiation  which  is  essential  for  the  proper  design  of  buildings, 
solar  energy  systems  and  a  good  evaluation  of  thermal  environment 
within  buildings. 

On  the  basis  of  the  models,  more  accurate  values  of  diffuse  radi¬ 
ation  can  be  obtained  which  will  be  helpful  for  the  prediction  of 
clearness  indices  I<T.  Clearness  index  is  an  important  factor  in  the 
design  of  solar  energy  systems.  A  greater  value  of  I<T  means  less 
diffuse  radiation  and  more  of  direct  radiation  favoring  concentrat¬ 
ing  solar  systems  while  non-concentrating  devices  are  needed  for 
locations  with  low  7<r.  Clearness  index  is  also  an  indicator  for  the 
installation  of  PV  panels-  the  optimal  tilt  angle  on  fine  days  is  large 
while  for  cloudy  days  the  tilt  angle  is  kept  low.  The  level  of  7CT  affects 
the  performance  of  solar  collectors  and  the  installation  of  PV  panels. 
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